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FUNDAMENTALS  OF  NON-NEWTONIAN  LIQUID  FLOW  ON  SURFACES 
1.  INTRODUCTION 

Chemical  agents  are  delivered  to  targets  in  a  number  of  different 
ways.  A  common  feature  of  the  delivery  protocols  is  that  the  agents 
originally  constitute  a  large  bulk  of  liquid,  which  is  then  somehow 
fragmented  or  atomized,  resulting  in  a  cloud  or  array  of  small  liquid 
drops.  In  order  to  ensure  that  the  drop  sizes  are  large  enough  (1mm  to  3mm) 
that  evaporation  is  not  dominant,  relatively  large  amounts  of  polymer 
are  added  to  the  bulk  liquid.  The  presence  of  the  polymer  directly 
affects  the  process  of  break-up  into  droplets,  and  thus  control  the  size 
distribution  of  the  drops. 

As  the  drops  are  propagated  through  the  environment  or  fall  under  the 
influence  of  gravity,  they  respond  to  change  in  the  ambient  temperature, 
pick  up  contaminants  from  the  atmosphere  and  lose  mass  due  to  evaporation. 
Eventually  they  strike  targets,  generally  solid  surfaces,  adhere  to  these 
surfaces,  begin  to  spread. 

In  order  to  devise  a  procedure  for  the  protection  from  and  clean-up 
of  these  agents,  one  must  understand  what  controls  the  spreading  process 
and  in  particular  the  parameters  that  govern  the  rate  of  spreading.  In 
effect,  the  clean-up  process  is  a  spreading  process  in  reverse;  the  removil 
of  a  drop  depends  on  similar  parameters  to  those  which  control  spreading. 

In  the  present  work  we  investigate  the  spreading  of  nou-Newtonlan 
drops  on  smooth  surfaces.  The  spreading  rate  thus  depends  on  the  surface 
tension  of  the  liquid-air  interface,  the  contact-angle  characteristics 
of  the  solid-liquid-air  system  and  the  material  properties  of  the  liquid. 

In  particular  we  assess  the  effects  on  spreading  of  viscoelastic  properties 
of  the  liquid  and  examine  the  influence  of  gravity. 


The  general  features  of  the  spreading-drop  model  are  described  in 
the  body  of  this  report,  together  with  a  presentation  and  discussion  of 
the  results  of  the  calculations.  The  details  of  the  calculations  are 
given  in  the  Appendix. 

2.  MECHANICS  OF  SPREADING 

Drops  that  strike  a  smooth  solid  spread  through  two  joint  mechanisms, 
contact-angle  "pull"  and  capillary  forces.  Both  contact-angle  "pull"' 
and  capillary  spread  are  included  in  the  present  model. 

2.1  Contact -Angle  ”Pull" 

Figure  1  shows  typical  measurements  of  contact  angle  6  versus  contact 

line  speed  u_  for  liquid-gas  systems.  When  liquid  displaces  gas  u„  >  0 
U*  CL 

and  angle  Increases  with  speed.  The  advancing  angle  6^  is  the  limiting 
angle  for  u—  0.  Thus,  a  drop  placed  on  a  solid  with  Initial  angle 
0O  >  9^  corresponds  to  a  contact  line  with  positive  (outward)  speed  and 
the  drop  is  pulled  outward  at  the  contact  line. 

2.2  Capillary  Spread 

If  the  drop  hits  the  solid  and  9q  ■  0^,  the  drop  will  spread  as  long 
as  the  drop  does  not  have  the  shape  of  a  stable  meniscus.  Surface  tension 
modifies  the  shape,  steepens  the  contact  angle  and  spreading  results. 

3.  RHEOLOGY 

There  are  three  general  features  common  to  the  motion  of  non-Newtonian 
fluids:  stress  relaxation,  normal  stresses  and  shear  thinning.  In  the 
present  formulation  we  employ  a  one  parameter  family  of  generalized  Maxwell 
fluids,  each  having  a  single  relaxation  time  x.  These  models  include  all 
three  non-Newtonian  effects  while  limiting  cases  suppress  either  stress 
relaxation  or  shear  thinning. 
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Figure  1 


Illustration  of  how  contact  line  speed  varies  with 
contact  angle 


3.1  Stress  Relaxation 


Stress  relaxation  is  the  time-dependent  adjustment  to  rapid  changes 
in  stress,  the  time  scale  (scales)  dependent  on  the  relaxation  time 
(times)  of  the  material. 

3.2  Normal  Stresses 

Normal  stresses  are  stresses  that  develop  in,  say,  pure  shear  flows 
in  directions  orthogonal  to  the  applied  force  and  are  generated  by  the 
nonlinear  behavior  of  the  material. 

3.3  Shear  Thinning 

The  nonlinear  behavior  further  allows  the  decrease  of  the  apparent 
shear  viscosity  p  with  shear  rate  y. 

4.  THE  MODEL  OP  THE  DROP 

The  drops  of  immediate  Interest  are  ones  that  have  small  initial 
angle  6^  and  spread  indefinitely  (to  final  angle  zero).  We  thus  use  an 
approximate  theory  valid  for  drops  whose  shapes  always  have  small  slopes, 
the  so-called  lubrication  approximation.  The  method  involves  a  systematic 
asymptotic  analysis  which  shows  that  inertial  effects  in  the  liquid  ate 
negligible.  More  importatly,  it  shows  that  if  shear  thinning  behavior 
is  present,  then  it  is  more  important  than  stress  relaxation  or  normal 
stresses.  Thus,  in  this  case  the  main  rheological  data  required  is  the 
zero  shear  rate  viscosity  and  the  changes  of  y  with  shear  rate.  On  the 
other  hand,  if  the  model  uses  a  constitutive  relation  for  which  shear 
thinning  (and  hence  normal  stresses)  are  absent,  the  stress  relaxation 
dominates. 

The  present  model  thus  incorporates  viscous,  rheological  and  surface 
tension  effects  in  the  spreading  process.  The  result  of  lubrication 
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analysis  is  a  simplified  set  of  nonlinear  partial  differential  equations 
(plus  initial  aud  boundary  conditions)  that  require  numerical  methods  for 
solution.  This  system  contains  six  non-dimensional  pa.  jeters. 


5.  PARAMETERS 

The  dimensional  quantities  that  enter  the  model  are  as  follows: 

P  density  of  the  liquid, 

zero-shear-rate  viscosity  of  the  liquid. 

p  surface  tensioh  of  the  liquid-gas  interface. 

K  reciprocal  of  the  slope  at  u^  ■  0  for  the 

characteristic  6  versus  u„  as  shown  in 

CL 

Figure  1,  k  has  units  of  velocity. 

0^  advancing  contact  angle  of  ths  drop. 

t  relaxation  time  of  the  liquid, 

g  magnitude  of  the  gravitational  acceleration. 

V  volume  of  the  drop, 

a^  initial  radius  of  the  drop. 

6q  initial  contact  angle  of  the  drop. 

The  above  ten  dimensional  quantities  form  into  six  non-dimensional  groups 
that  characterize  the  spreading  of  drops.  These  are  us  follows: 


capillary  number 


relaxation  parameter  #1 


relaxation  parameter  #2 


Bond  number 


B 

e. 


2 

Pga0 


eA/e0 


final  contact  angle 


V  - 


»oe0 


volume  of  the  (axisymmetric)  drop 


‘o8o 


volume  per  unit  depth  of  the  (two-dimensional) 
drop. 

The  capillary  number  expresses  the  relative  importance  of  contact- 

angle  "pull",  measured  by  speed  k9q,  to  capillary  spread,  measured  by 

06q/v<q.  The  relaxation  parameter  #1  measures  the  relaxation  time  t  versus 

capillary  effects.  The  relaxation  parameter  #2  measures  the  relaxation 

time  t  for  shear  thinning  versus  capillary  effects.  Notice  that  A/e  -  0Q 

and  since  6q  is  very  small,  stress  relaxation  effects  are  small  compared 

to  those  due  to  shear  thinning  except  when  the  latter  are  absent.  The 

Bond  number  measures  the  effects  of  gravity  (compared  to  surface  tension) 

when  the  plate  is  horizontal.  We  allow  the  possibility  that  the  drop 

ceases  spreading  at  final  angle  0p«  If  6^  -  0,  the  drop  spreads  indefinitely. 

In  order  to  see  the  ranges  of  values  taken  by  the  various  parameters, 

we  write  down  "typical"  values  for  the  dimensional  quantities.  We  use 

o  *  -4  3 

a^  -  0.25  cm,  8Q  »  0.1  rad  *  5.7  ,  V  ■  1.5  x  10  cm  ,  a  ■  30  dynes/cm, 

3 

yQ  -  5  poise  and  p  ■  l  gm/cm  •  We  have  no  Information  on  the  contact  angle 
constant  k;  we  take  tc  ■  0.1  cm/ sec.  The  relaxation  time  t  can  be  obtained 
from  data  on  apparent  viscosity  versus  shear  rate  and  equation  (2.44)  of 
the  Appendix  if  we  rewrite  it  in  dimensional  form.  The  apparent  viscosity 
li  of  our  model  takes  the  form 
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i  j.  2"2 

1  +  T  Y 


(1) 


where  is  the  zero-shear-rate  viscosity  and  Y  is  the  shear  rate.  From 
the  CSL  data  given  to  us,  we  estimate  for  Diethylmalonate  thickened  with 
2.1%  Polymethylmethacrylate  that  t  ■  1.0  sec;  for  Diethylmalonate 
thickened  with  5.27.  Copolymer,  t  *  0.013  sec;  for  Diethylmalonate  thickened 
with  9.5%  Elvacite  2041,  t  =  0.0063  sec;  for  Methyl  Salicylate  thickened 
with  4.5%  Copolymer,  t  -  0.010  sec.  *  us,  a  practical  upper  limit  for  t  is 
about  1  sec.  With  this  estimate  we  find  that  "typical"  values  for  the 
parameters  are  as  follows: 


C  «  1.6 

X  «  2.5  x  10“2 

e  «  2.5  x  10 (2) 

B  «  2.1 

V  *  1.0 

Clearly,  the  estimates  (2)  are  crude  and  each  parameter  haB  a  range  of 
possible  values,  we  shall  use  e  «  5.0,  a  very  large  value  for  the  calcu¬ 
lations  in  order  to  exaggerate  the  effects  of  viscoelasticity. 

6.  RESULTS 

We  use  the  spreading  of  a  Newtonian  liquid  drop  as  the  standard 
against  which  to  compare  results  for  non-Newtonian  cases.  We  thus  first 
present  resulta  for  this  case.  Since  spreading  rates  depend  explicitly 
on  time  we  shall  define  a  non-dimensional  time  t  in  terms  of  the  dimensional 

A 

time  (sec)  t  as  follows: 


In  terms  of  the  typical  values  of  parameters  given  in  Section  5,  t  ~  1 

ic 

corresponds  to  physical  time  t  about  2.5  seconds, 

6,1  Two-Dimensional  Newtonian  Drops 

The  discussion  here  relates  to  a  Newtonian  constitutive  model  in  which 
the  viscosity  is  constant,  but  y0  should  be  regarded  for  purposes  of 
comparison  as  the  zero-shear-rate  viscosity  of  the  real  liquid. 

Figures  2a  and  2b  show  the  shape  of  the  drop  at  various  Instants  of 
time  t  for  the  case  V  -  0.75,  C  -  1,  X  -  e  -  B  -  0.  In  Figure  2a  the 
final  angle  0F  -  0.5  and  so  the  drop  shape  approaches  a  static  equilibrium 
as  t  -*■  In  Figure  2b,  the  final  angle  0p  ■  0  and  so  the  drop  spreads 
indefinitely.  These  Figures  refer  to  the  early  stages  of  the  drop's  spread. 
By  contrast  we  illustrate  in  Figures  3a  and  3b  the  shape  of  the  drop  when 
substantial  spread  has  occurred.  In  these  two  Figures  we  have  taken 
0p  ■  0,  C  ■  1,  X  ■  e  *  B  “  0,  and  we  show  the  shape  of  the  drop  when  t  ■  0 
and  when  t  -  t*2\  where  t*2^  is  the  time  for  the  drop  to  spread  to  a 
size  twice  its  initial  radius.  We  shall  call  a(t)  the  radius  of  the  drop 
even  though  in  two  dimensions  it  is  really  the  half-width  of  the  drop. 

Note  that  a(0)/aQ  -  1  by  definition.  In  Figure  3a  we  have  taken  V  -  0.75, 
while  V  -  1.0  in  Figure  3b. 

Graphs  such  as  the  above  give  qualitative  information  but  quantitative 
predictions  must  be  extracted. 

First,  we  examine  cases,  6p  «  0,  for  which  the  drop  spreads  indefinitely. 
We  set  X  ■  c  -  B  -  0  and  examine  t^  as  functions  of  both  C  and  V. 

(O') 

Figure  4  shows  '  versus  V  for  several  values  of  C.  The  spreading  rate 
depends  strongly  on  the  volume  of  the  drop  (for  fixed  initial  contact 
angle  6q),  large  drops  spread  faster  than  smaller  drops.  Also,  for  a 
fixed  volume  of  liquid  the  time  it  takes  to  double  the  initial  radius 
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Figure  3a:  Tvo-dimanslonal  drop,  shapes  at  t  *  0  and  t  -  t<2)  for 
V  -  0,7  j,  C-l,  e-X-B-0 
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of  the  drop  is  larger  the  smaller  the  value  of  C.  C  decreases  with  surface 
tension  a  and  Increases  with  zero- shear -rate'  viscosity. 

Figure  5  shows  the  doubling  time  tv  as  a  function  of  the  capillary 
number  C  in  the  case  V  ■  1.  Figure  6  gives  the  radius  a(t)  versus  t  for 
various  values  of  V,  here  for  C  ■  1.  Figure  7  shows  the  radius  a(t)  versus 
t  for  various  values  of  C,  here  for  V  *  1*  Thus,  Figures  4  and  5  summarize 
the  data  of  Figures  6  and  7. 

In  Figure  8  we  extend,  on  a  different  scale,  the  variation  of  a(t) 
with  t  in  the  case  V  ■  1,  C  ■  1. 

We  now  examine  the  case  for  ^  0.  Although  the  case  0?  -  0  is  of 
most  interest  for  CSL,  we  used  cases  for  i  0  as  tests  of  the  numerics 
and  we  present  some  of  the  results,  here* 

We  confine  our  attention  to  1  ■  e  ■  B  ■  0  and  take  -  0.5.  Figure  9 
shows  the  radius  of  the  drop  as  a  function  of  time  for  various  V  with 
C  ■  1.  The  vertical  lines  are  the  asymptotes  valid  for  t  ®.  There  seems 
to  be  little  qualitative  difference  between  cases  e?  -  0  and  J*  0  until 
one  approaches  times  large  enough  that  the  finite  equilibrium  is  closely 
approached. 


6.2  Axisymaatric  Newtonian  Drops 

Figures  10-14  give  for  the  axisytnetr-ir  drop  the  information  equivalent 
to  Figures  4-8  of  the  two-dioanalonal  drop.  The  same  trends  persist  here 


as  in  that  case*  Spreading  is  more  rapid  for  large  drops  that  have  large 
capillary  numbers.  However,  we  do  see  that  axisymnetric  drops  spread 
more  slowly  than  their  two-dimensional  counterparts.  For  example  if 
C  *  V  “  1,  t^  ■  13.0  for  the  two-dimensional  drop  while  -  27.2 

for  the  axisyntnetric  drdp. 
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Figure  S:  Two-dimensional  drop,  variation  of  doubling  time  with  C 


for  various  B,  V  »  1,  e  “  X  ■  0 


1.0 


1.5 


2.0 


Figure  6:  Two-dimensional  drope  variation  of  radius  with  t  for  various 

V,  with  C  “  1,  9  «•  0,  e  ■  X  -  B  ■  0 

F 
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rate  of  change  of  radius  for  various 

■  0 


6.3  Two-Dimensional  Viscoelastic  Drops 


The  constitutive  model  used  here  contains  a  single  relaxation  time 
which  enters  the  governing  equations  through  two  parameters,  the 
relaxation  parameters  X  and  e.  Recall  that  X  is  formally  much  smaller 
than  e  by  a  factor  6q. 

We  now  present  spreading  curves  for  non-Newtonian  cases  and  compare 
these  with  the  corresponding  Newtonian  cases.  Figure  15  shows  the  radius 
of  the  drop  as  a  function  of  t  for  three  cases  with  C  ■  V  ■  1,  B  ■  0, 

6^  -  0,5,  There  is  the  Newtonian  case  e  ■  X  ■  0,  a  stress-relaxation 
case  e  ■  0,  X  «  0.5  and  a  shear  thinning  case  e  »  5.0,  X  -  0,  Firstly, 
we  see  that  for  equal  values  for  the  material  relaxation  time  t  and  for, 
say,  ■  0,1  that  £  -  10X  and  the  effect  of  shear  thinning  is  much  larger 
than  that  of  stress  relaxation.  Shear  thinning  causes  the  drop  to  spread 
faster  than  the  corresponding  Newtonian  drop  given  that  the  two  have 
equal  zero -shear-rate  viscosities.  For  example  in  the  present  example 
the  viscoelastic  drop  reaches  150%  of  its  Initial  radius  when  t  ■  4.0 
while  the  corresponding  drop  for  X  ••  0  reaches  this  position  at  t  *  6.4, 

If  stress  relaxation  only  is  considered,  the  drop  spreads  more  slowly  and 
attains  150%  of  its  initial  radius  ut  t  "  6,8.  Rather  than  pursue  ranges 
of  parameters  here,  we  turn  to  the  axisymmetric  drop  and  confine  our 
attention  to  the  principal  case,  X  «  0,  e  0,  in  which  shear-thinning 
only  is  present. 

6.4  Axisynmetric  Viscoelastic  Drops 

Figures  16-18  give  results  for  drops  that  have  significant  shear¬ 
thinning  effects  e  «  5.0,  vanishingly  small  stress  relaxation  X  =  0  and 
final  contact  angle  0  «  0  so  spreading  continues  indefinitely.  Gravity 
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Figure  16;  Axisynmetric  drop,  variation  of  radius  with  t  for  shear-thinning 


is  i(iored,  B  *  0,  and  a  range  of  capillary  numbers  C  and  drop  volumes 

V  is  examined. 


Figure  16  gives  the  radius  a(t)  of  the  drop  as  a  function  of  time 
for  B  ■  A  *  0  and  G  »  V  “1.  There  are  two  cases  shown,  one  for  the 

Newtonian  drop  with  e  -  0  and  one  for  the  viscoelastic  drops  with  e  ■  5,0. 

We  see  that  t^^  decreases  by  about  157.  due  to  shear  thinning. 

We  can  compare  the  viscoelastic  cases  to  the  Newtonian  cases  by 
(21 

defining  a  ratio  R'  , 

t<2> 

(2)  VISCOELASTIC  . 

R  -  (2)  .  W 

'"NEWTONIAN 


The  quantity  R^1*5^  is  the  similar  ratio  for  a  50%  increase  in  radius. 

Figure  17  shows  R^  and  R^1,5^  for  the  case  B  -  0,  V-C-l  where 
the  right  aide  of  the  figure  shows  shear-thinning  effects  only  e  i1  0, 

X  -  0  and  the  left  side  shows  stress  relaxation  effects  only,  c  -  0,  X  +  0. 
On  one  hand  we  again  see  that  shear  thinning  has  a  larger  effect  than  does 
stress  relaxation.  On  the  other  hand  there  is  a  "catch-up"  phenomenon 
present.  The  largest  effect  on  spreading  of  viscoelasticity  occurs  at 
early  times  t.  However,  not  only  does  viscoelasticity  have  weaker  effects 
later  but  the  differences  in  spread  rates  become  smaller  as  time  passes. 
Notice  that  R^  is  nearer  to  unity  than  is  R^’^.  Thus,  viscoelastic 
drops  that  have  had  a  long  time  to  spread  have  attained  nearly  the  same 
radii  as  their  Newtonian  counterparts.  Figure  18  shows  equivalent  informa¬ 
tion  for  the  case  B  -  0,  C  “  0.1,  V  -  1.  Notice  that  the  effects  of  visco¬ 
elasticity  are  much  smaller  here  for  C  -  0.1  than  for  the  earlier  case  for 
G  -  1.0. 
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When  a  drop  spreads  on  a  smooth  horizontal  plane,  gravity  modifies 


the  shape  and  spreading  characteristics  of  the  drop. 

Figure  19  shows  a  two-dimensional  Newtonian  case  for  C  *■  1,  V  -  0.75 
and  B  ■  10  at  various  times  t.  A  comparison  of  these  curves  with  those 
of  Figure  2b  for  the  equivalent  case  with  B  -  0  shows  that  gravity 
steepens  the  interface  near  the  contact  line  which  leads  to  more  rapid 
spreading.  Similar  behavior  can  be  observed  in  Figures  20a  and  20b,  which 
are  the  direct  analogs  of  Figures  3a  and  3b  respectively  for  the  case  B  ■  10. 

Figures  21  and  22,  respectively,  show  a(t)  versus  (t)  for  various 
values  of  B  for  two-dimensional  and  axisymmetrlc  spread.  Notice  for 
B  ■  1.0  that  ■  0.76  for  the  two-dimensional  drop  and  ■  0.84  for 

the  ax 1 symmetric  drop  so  that  there  are  comparable  accelerations  of 
spreading  in  both  geometries  due  to  gravity.  Effects  of  gravity  are  also 
Indicated  in  Figures  8  and  14,  which  give  long-time  spreading  behavior, 
and  in  Figures  5  and  11,  which  show  how  doubling-time  varies  with  capillary 
number • 

6.6  Effects  of  Gravity  on  Newtonian  Drops  on  Tilted  Surfaces 

A  drop  sits  on  a  plane  tilted  an  angle  3  to  the  horizontal.  The 
component  g  cos  8  of  gravity  normal  to  the  plate  leads  to  a  Bond  number 
B  cos  @  whose  effect  was  discussed  in  Subsection  6.5.  The  component  gain  3 
of  gravity  "slides"  the  drop  bodily  down  the  plate  as  discussed  by  Hocking* 
for  the  two-dimensional  case  with  C-*  . 

The  wetting  properties  for  C  <  »  shown  in  Figure  1  allow  the  drop  to 
make  an  initial  adjustment  in  its  shape  to  accommodate  to  a  steady 

*k 
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Figure  20b:  Two-dimeruional  drop,  shapes  at  t  »  0  and  t  ■  t 
V  -  l,  C  -  1,  B  -  10,  e  -  X  -  0 


for 


(2) 


t 
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translation*  Tills  accommodation  takes  place  as  follows.  Initially,  the  rear 
contact  angle  is  ancj  the  forward  contact  angle  is  q^r®NT.  If 

0*  <  ©a  (advancing)  angle)  and  9°  <  0^  (receding  angle)  then  both 

contact,  lines  move.  If,  according  tc  Figure  1,  the  corresponding  contact¬ 
line  speeds  are  equal,  then  the  drop  "slides"  down  the  plane  with  constant 
shape.  If,  for  example,  .the  front  contact  line  initially  moves  faster 
than  the  back  one,  the  hydrodynamics  causes  the  front  angle  to  become  more 
shallow  slowing  that  contact  line  and  the  back  angle  to  become  more  shallow 


speeding  that  contact  line.  The  result  tends  toward  the  steady  sliding 
case.  Clearly,  other  cases  also  "rearrange"  and  approach  the  steady 
sliding  case  as  along  as  0R  >  0. 

The  net  result  of  these  observations  Is  that  the  component  of  gravity 
along  the  plate  has  only  indirect  effects  on  spreading  characteristics 
since  it  merely  adjusts  the  contact  angles.  The  component  of  gravity 
normal  to  the  plate  directly  affects  spreading  as  discussed  in  Subsection  6.5. 


7.  CONCU  >IONS 

The  characteristics  of  Newtonian  liquid  drops  spreading  on  smooth 
horizontal  solid  planes  can  be  characterized  by  the  non-dimensional  groups 
V  and  C  listed  in  Section  5.  The  intrinsic  spreading  and  capillary  effects 
are  contained  in  the  capillary  number  Cj  C  is  large  if  the  zero-shear-rate 
viscosity  is  large  or  the  surface  tension  is  small.  The  bulk  features 
of  the  drop  are  contained  in  the  volume  V.  The  drop  spreads  more  rapidly 
if  C  and  V  are  increased. 

The  viscoelastic  effects  of  the  drop  are  characterized  by  the.  parameters 
1  and  e  listed  in  Section  5.  Both  X  and  e  measure  the  (single)  relaxation 
time  of  the  material  though  these  are  interpreted  as  stress  relaxation 
and  shear-thinning  measures,  respectively.  For  small  contact  angles 
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e  »  \  so  that  shear  thinning „  if  present:,  is  much  more  important  than 

stress  relaxation. 


We  find  that  shear  thinning  is  the  dominant  viscoelastic  effect  in 
thin  drops.  It  has  its  greatest  effects  at  small  times  (the  drop  begins 
spreading  at  t  ■  0) .  For  example  shear  thinning  can  accelerate  the  spreading 
rate  by  50%  after  a  few  seconds  if  e  ■  5.0.  However,  there  is  a  "catch-up" 
phenomenon  that  makes  the  spreading  rate  for  long  times  very  nearly  that 
appropriate  to  the  Newtonian  drop  having  the  same  zero-shear-rate  viscosity. 

The  main  conclusion  one  draws  from  the  present  study  is  that  the 
principal  rheo log leal  property  of  interest  in  the  spreading  mechanics  of 
drops  on  smooth  surfaces  is  Pq,  the  zero-shear-rate  viscosity.  Thus  for 
large  ranges  of  material  of  Interest  to  CSL,  polymer  additives  greatly  modify 
Pq  and  spreading  is  controlled  by  pg. 

The  second  most  important  rheological  property  is  the  variation  p(Y) 
of  shear  viscosity  with  shear  rate  Y.  In  the  calculations  presented  we 
used  the  estimate  t  ■  1  sec  for  the  maximum  value  of  material  relaxation 
time.  Under  most  practical  circumstances  for  CSL  e^O.Ol  though  it  can 
in  extreme  cases  get  as  large  as  5.0.  Even  for  b  -  5.0  we  found  rather 
modest  acceleration  in  the  spreading  rate  (compared  to  the  case  e  *  0  with 
the  same  Pq) ;  small  values  of  e  yield  only  minor  accelerations.  Hence, 
viscoelasticity  affects  spreading  in  a  calculable  way  but  its  magnitude  is 
seemingly  small. 

The  effect  of  a  vertical  gravity  on  the  spreading  of  a  drop  on  a  smooth 
horizontal  plate  is  characterized  by  the  Bond  number  B  defined  in  Section  5. 
Vertical  gravity  distorts  the  shape  of  the  liquid-air  interface,  steepens 
the  contact  angle  and  thus  accelerates  the  spreading  process.  The  effect 
is  appreciable  for  large  drops  and  negligible  for  small  drops. 


A  drop  on  an  inclined  plane  tilted  an  angle  6  to  the  horizontal  has 
a  component  of  gravity  normal  to  the  plate  that  affects  spreading,  as  above, 
through  a  Bond  number  BcosS.  The  component  of  gravity  down  the  plate  will 
result  in  the  drop  having  steeper  contact  angles  in  front  and  shallower 
angles  in  the  rear.  The  relations  between  these  changes  affect  the  spreading 
somewhat  but  principally  lead  to  the  drop  moving  bodily  down  the  plate  but 
not  experiencing  major  changes  in  spreading  rates. 
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APPENDIX 
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1.  Formulation 

2.  Lubrication  Approximation 

3.  Evolution  Equation 


4,  Solution  Procedure 


1.  Formulation 


We  consider  a  drop  of  viscous  liquid  on  a  smooth  rigid  horixontal 

^  ^  ^ 

plane.  We  use  a  Cartesian  coordinate  system  (x  ,y  ,z  )  for  a  two-dimensional 

*it  fc  * 

drop,  and  a  cylindrical  polar  coordinate  system  (r  ,  ,z  )  for  an  axi- 
synmetrlc  drop.  In  both  cases  we  take  the  rigid  plane  to  coincide  with 

‘M  1 

the  plane  z  ■  0,  with  the  z  axis  pointing  vertically  upwards. 

In  the  case  of  a  two-dimensional  drop  we  consider  the  motion  to  be  in  the 
x  ,z  -plane,  with  all  quantities  independent  of  y  .  Equations  and  rela¬ 
tions  pertaining  to  the  two-dimensional  case  are  designated  by  the  letter  P. 
For  the  axisymnetric  drop  all  quantities  are  independent  of  the  azimuthal 

*k 

coordinate  <j>  ;  this  case  will  be  designated  by  the  letter  A. 
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with 


W  - 0  ; 

the  contact  angle  0^  i3  given  by 


(1.2A) 


dh 


* 


tan  9q  ■  ■ 

and  the  volume  is 

a. 


dr 


£  <v  . 


(1.3A) 


V  m  2tt  I  r  hQ(r  jdr  . 


<1.4A) 


Drop  Dynamics 

At  time  t  >  0  we  denote  by  a  (t  )  the  half -width  of  the  drop  in  the 
two-dimensional  case  and  its  radius  in  the  axisymmetric  cose*  Thus  we 
have  the  initial  condition 


a  (0)  «  a 


0 


We  take  tha  shape  of  the  drop  at  time  t  to  have  the  form 
ft  "  h  (x  ,fc  )  ,  -a  (t  )  <  x  <  a  (t  ) 


with 


*  *  *  *  *  *  *  * 
h  (-a  (t  ),t  )  ••  h  <&  (t  ),t  )  -  0 


in  the  two-dimensional  case,  and 

*,*.#*.  „  *  *  * 
z  -  h  (v  ,t  )  ,  0,  <_r  <  a  (t  ) 


(1.5) 

(1.6P) 

(1.7P) 

(1.6A) 


with 


h  (a  (t  ),t  )  -  u 


(1.7A) 


in  the  axisymnetric  case.  The  contact  angle  at  time  t  is  denoted  by 

^  ^  yt. 

6  *  6  (t  ),  and  is  given  by 


Appendix 


49 


v,*.  .«.«  .V* 


wl'«. 


(1.8P) 


*  3hv  ,  *  *.  *v  3h*  ,  *K  *V 

tan  e  -  —  (-a  (t  ),t  )  - - -  (a  (t  ),t  ) 

ax  a  x 


„*  3h  ,  *,  *  *y 
tan  0  -  -  — r  (a  (t  )  ,t  ) 

3r 


(1.8A) 


respectively.  Conservation  of  the  volume  of  the  drop  over  time  gives  the 


additional  constraint 

*.  . 

a  (t) 


*  f  .  *  *  *  .  * 
V  «  I  h  (x  ,t  )dx 

J  *  .  *. 

-a  <t  ) 


(1.9P) 


V*  ■  2ir 


*  * 
*  <t  ) 


/*  *  *  *  * 

r  h  (r  ,t  )dr  . 


(1.9A) 


iKr  ^ 

The  rate  of  change  of  the  quantity  a  at  any  instant  t  is  taken,  on 
empirical  grounds,  to  be  proportional  to  the  difference  between  the  angle 
9  and  the  advancing  contact  angle  6^.  The  latter  is  the  static  equilibrium 
angle  when  the  drop  is  on  the  point  of  spreading.  The  implied  relation  has 


the  form 


K(0*-  0  ) 
dt  A 

where  k  >  0  is  an  empirically  determined  constant,  and  where  9^  >.  0. 


(1.10) 


Equations  and  Boundary  Conditions 

The  motion  is  governed  by  the  Navier-Stokes  and  continuity  equations, 


/8s  *  *  *\ 


it  it  it  it  A 

-  V  p  +  V  *S  -  pgz 


V*‘v*  -  0 


(1.11) 


(1.12) 


' ft  ^ 

where  p  is  the  density,  v  *  (u  ,0,w  )  is  the  velocity  vector,  p  is  the 

it 

pressure,  S  Is  the  extra-stress  tensor,  g  is  acceleration  due  to  gravity, 
Appendix  50 

-v  v.  '.v  c.  w*  ■  -*■  ■ 


und  z_  is  unit  vector  in  the  upward  vertical  direction*  A  constitutive 
relation  between  the  stress  and  the  deformation— rate  will  be  given 
below. 

The  boundary  conditions  are: 

(i)  The  normal  velocity  component  is  zero  on  the  rigid  plane i 

w*  «  0  on  z*  “  0  ;  (1.13) 

(ii)  As  explained  by  Dussan  V.  and  Davis  [1],  the  usual  no-slip  condition  at 
the  rigid  boundary  needs  to  be  modified  to  avoid  the  appearance  of  a 
singularity  at  the  contact  line*  Following  Greenspan  [2]  we  take  the 
condition  to  be 

(Y  Z/n0>Si3  -  h  u  on  z  «  0  (1.14) 

where  PqIs  the  zero -shear-rate  viscosity,  Y  is  a  slip  length  (distance 

‘k 

from  the  contact  line  over  which  slip  takes  place)  and  S^3  represents  the 

k  k 

x  z  -component  of  extra  stress  in  the  two-dimensional  case,  and  the 
*  *  .  , 

r  z  -component  in  the  axisymmetric  case. 


(Hi)  The  kinematic  condition  at  the  free  surface  is 


*  3h*  ,  *  3h*  *  .  *>  *  *. 

w  »  — -y  +  u  — -  on  z  ■  h  (x  ,  t  ) 

3t  3x 


(1.15P) 


or 


ju  Ju 

*  3h  .  *  3h 

w  -  -*  +  u 

3t  8r 


on  z  ■  h  (r  ,t  ) 


(1.13A) 


(iv)  The  dynamic  boundary  condition  at  the  free  surface  is 
**  *  ft* 

“  tP  Je  +  S  *S.  “  o  n  on  z  “  h  (1.16) 

where  [p  ]  denotes  the  pressure  difference  across  the  interface,  a  is  the 
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surface  tension,  n  is  the  outward  unit  normal  to  the  surface,  and  2H 
is  the  mean  curvature  of  the  interface. 

Rheology 

A  constitutive  equation  that  is  suitable  for  polymer  solutions  is  the 
generalized  Maxwell  model,  ' 

* 

*  *  l.w* 

2  +  T  \  — y  +  (V  -V  )S  +  r-  <«  -g  -  S  •“  ) 

fat 

i  *  •*  |  •*  ..  . 

-  2  e<2  *X  +  X  2  >  f  -  vl  (1.17) 


where 


.•*  *  *  „  *  *  T 

l  -  ?  v  +  (V  v  ) 


*  *  *  ,  *  *  T 

Uj  -  V  V  -  (V  v  ) 


(1.18) 


(1.19) 


t  is  the  relaxation  time  of  the  liquid  and  0  is  a  number,  which  in  practical 
situations  can  range  between  -1  and  +1,  (See,  for  example,  Petrie  [3].) 
Rather  than  analyze  the  whole  range  of  values  of  8,  we  shall  consider  only 
two  special  cases,  namely  0  ■  0  and  0  ■>  1.  When  0  -  0  the  model  (1,17) 
reduces  to  the  well-known  corotatlonal  Maxwell  model.  In  steady  uni¬ 
directional  shear  flow  this  model  yields  shear  thinning  and  both  first 
and  second  normal  stress  differences,  and  yields  stress  relaxation  in 
unsteady  simple  shear.  When  0  «  1,  (1,17)  reduces  to  the  upper  convected 
Maxwell  model.  Here  there  is  no  shear  thinning  in  simple  shear,  the  vis¬ 
cosity  remaining  constant  at  its  zero-shear- rate  value,  but  first  and 
second  normal  stress  differences,  as  well  as  stress  relaxation,  are 
present.  By  studying  these  two  models  ws  have  the  possibility  of  assessing 
the  various  rheological  effects  on  spreading;  In  particular  the  differences 
between  the  0  ■  0  and  0  «  1  case  will  distinguish  the  influence  of  shear 
thinning. 
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2.  Lubrication  Approximation 

We  proceed  on  the  basis  of  the  assumption  that  the  initial  angle  9^ 
is  very  small.  This  enables  us  to  use  the  lubrication  approxi¬ 
mation,  in  which  all  quantities  are  appropriately  scaled,  and  then  the 
equations  und  boundary  conditions  are  expanded  in  powers  of  0q.  The  firBt- 
order  problem  is  what  is  retained  in  the  limit  0^  -+•  0. 

A  dimensionless  time  t  is  defined  by 

t-^^Vo  '  <2*l> 


and  dimensionless  coordinates  (x,z)  in  the  two-dimensional  problem,  and 
(r,z)  in  the  axisymmetric  problem,  are  defined  by 

x  -  x*/&Q  ,  t,  -  z*/(aQ0o)  ,  r  -  r*/aQ  .  (2.2) 

The  dimensionless  shape  of  the  drop  becomes  h(x,t)  or  h(r,z)  respectively, 
where 

h  "  h*/ (aoecP  *  (2.3) 

The  half-width  (radius)  of  the  drop  is  a(t),  given  by 

a(t)  -  a*(t*;/aQ  (2.4) 

and  the  contact  angle  is 

0(t)  -  0*(t*)/Oo  (2.5) 

with  the  final  equilibrium  (advancing)  contact  angle  given  by 

\  -  Veo  •  <2-6> 

We  also  define  a  dimensionless  volume  by 

V  -  V*  / ao 90  (2.7P) 

or 

v  =  ,v*/a  .  (2.7A) 


Appendix 


53 


Dimensionless  velocity  components  (u,w)  and  pressure  p  are  given  by 


U  -  U*/(K0O)  ,  w  =  w*/(ice2) 


P  <‘aO0O/wOlc) 


(2.8) 


Finally  we  have  a  dimensionless  stress  tensor  S  defined  by 

2  "  2*(a0/vi0K)  *  <2-9> 
Initial  state 

Nose  that  our  non-dimensionalization  gives  that  the  drop  has  unit 
half -width  (radios)  initially,  and  the  initial  contact  angle  3(0)  is 
also  unity. 

In  the  two-dimensional  case  the  initial  shape  has  the  form 

z  ■  hQ(x)  ,  -1  <  x  <  1  (2.10P) 


with  the  end  conditions 

ho<-D  -  VD  “  0 


(2. HP) 


Using  the  lubrication  approximation  6q  +0  we  find  the  condition  for  the 
initial  contact  angle  to  be 


(-1) 


dhc 

dx 


<-D  -  1 


(2.12P) 


hQ(l)  -  0  , 


(2.11A) 


at  the  contact  line 


<l)  -  -1  , 


(2.12A) 


-  hfr 


rh(r)dr 


(2.13A) 


The  shape  of  the  drop  at  time  t  in  the  two-dimensional  case  is 


z  ■  h(x,t)  ,  -a(t)  <  x  <  a(t) 


(2.1AP) 


h(-a(t),t)  -  h(a(t),t)  -  0 


(2.15P) 


while  in  the  axi 


Lc  case  it  1b 


z  »  h(r,t)  ,  0  <  r  <  a(t) 


(2.14A) 


h(a(t),t)  ■  0 


(2.15A) 


In  the  lubrication  approximation  the  contact  angle  0(t)  is  given  by 


0<t)  -  hx<-a(t),t)  hx(a(tj,t) 


(2.16P) 


0<t)  “  -  hr(a(t),t)  . 


(2.16A) 
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The  volume  conservation  condition  (1.9)  becomes 


or 


a(t) 

V  "  i  h(x,t)dx 
-a(t) 


a(t:) 

V  -2  TrJ  rh(r,t)dr  . 

0 


(2.17P) 


(2.17A) 


Equations  (1.10)  and  (2.16)  combine  to  give  the  following  differential 
equation  for  the  dynamics  of  the  drop: 


or 


"  C[~hx(a(t),t)  -  ep] 


-  Ct-hr(a(t),t)  -flp] 


where  C  1b  the  capillary  defined  by 

C  - 


V 


ere! 


(2.18P) 


(2.18A) 


(2.19) 


Equation  (2.18)  is  subject  to  the  initial  conditions 


a(0)  «  1 


(2.20) 


Equations  and  Boundary  Conditions 

In  the  lubrication  limit  0^  -*•  0  the  Navier-Stokes  and  continuity  equa¬ 
tions  in  the  two-dimensional  case  reduce  to 


-  "x  +  S13,x  ■  0 

(2.21P) 

-  p  -  B/C  -  0 

Z 

(2.22P) 

u  +  w  ■  0 

X  z 

(2.23P) 

where  Sjj  is  the  x,z  component  of  the  stress  tensor, 

Bond  number  defined  by 

and  where  B  is  the 
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B 


(2.24) 


pga2 

_ 0 

a 

Note  that  S^»  which  is  associated  with  shear  thinning  in  viscoelastic 
materials,  Is  the  only  stress  component  that  remains  in  the  reduced  system 
(2.2l)-(2,23).  Components  such  as  S^,  822*  associated, with  normal 

stress  differences,  are  not  explicitly  present  in  the  equation  of  motion 
in  the  lubrication  limit. 

The  corresponding  equations  in  the  axisymmetric  case  are 


-  p  +  S, »  "0 

*r  13, z 

(2.21A) 

-  p  -  B/C  -  0 

Z 

(2.22A) 

(ru)r  +  <rw)z  -  0 

(2.23A) 

The  boundary  conditions  are: 

(i) 

zero  normal  velocity  at  the  rigid  plane 

w  "  0  on  z  -  0 

(2.25) 

(ii) 

modified  slip  condition  at  the  rigid  plane 

2 

Y  Sj  j  «  hu  on  z  ■  0 

(2.26) 

where 

Y2  -  (Y*/a0V2  * 

(2.27) 

(iii) 

kinematic  interfacial  condition 

h  +  C(uh  -w)  *  0  on  z  -  h(x,t) 
c  X 

(2.28P) 

or 

h  +  C(uh  -w)  ■  0  on  z  -  h(r,t) 

C  IT 

(2.28A) 
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TS3KH.XJ..K.I  Tli.VS.Vv  x:  -T.' 


I 


ms*.**  ’r?  tv  *4;  •?  7.w> 


(iv)  dynamic  boundary  conditions  on  the  interface  are 
h  +  Cp  ■  0  on  z  ■  h(x,t) 

XX 


(2.29P) 


or 


and 


li  +  —  h  +  Gp  «  0  on  z  »  h(r , t) 
rr  r  r 


■  0  on  z  ■  h 


(2.29A) 


(2.30) 


Rheology 

Using  the  scalings  Indicated  above  we  find  that  the  constitutive 
relation  (1.17)  becomes,  in  dimensionless  form, 

|  +  \  ||  +  ecj  (v.V)g  +  \  <U|*g-B*f||)  -  |  B(g*Y  +  Y«g)J  -  l  (2.31) 
where  In  the  lubrication  approximation 


1  0  u. 

t  -  (  * 

u  0 


0  «u 


u  0 
z 


and  where  A,  £  are  relaxation  parameters  defined  by 
,3 


TO0 


0 


Vo  ’ 


(2.32) 


(2.33) 


8 


I 

i 

1 

w 

vf. 

ft 

1 


e  - 


T06o 

Vo 


Note  that 

A/e  -  0, 


(2.34) 


(2.35) 


The  form  of  (2.31)  shows  that  A  is  a  measure  of  stress  relaxation  while  e 
is  a  measure  of  shear  thinning.  Hence  (2.35)  implies  that  stress  relaxation 
only  becomes  Important  when  shear  thinning  is  absent. 
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k'j  -".  *i_j  -\  -r;  -■. 


i;h,;.\^v:\v\w;,i,.',',;w.uv'>1 


The  components  of  equation  (2.31)  are 
3S 


all  (  \ 

*11  +  X  —  +  60  |  ~uz*13  ‘  6Vl3}  -  0 

+  i  ec { u2(Su  - s33)  -  su«(su  +  s33)} 

SS33  f  1 

S33  +  X— +E°{Vl3-8u.S13>}-°  • 


3S13  1 

S13  +  X  It 


u 


(2.36) 

(2.37) 

(2.38) 


When  P  ■  1  (upper  convected  Maxwell  model)  we  obtain  ■  0  from  (2.38) 
and  then  (2.37)  reduces  to 

SS13 


%3  +  X  tc"  "  S 


<2.39) 

In  this  model  only  stress  relaxation  is  relevant,  and  there  is  no  shear 
thinning. 

When  (3-0  (corotational  Maxwell  model)  we  neglect  X  by  virtue  of  (2.33) 
and  the  system  (2.36) -(2.38)  reduces  to 


S11  -  eCuzSl3  "  0 


S13  +  2  ^z^ll-V  "  uz 


S33  +  eCuzS13  "  0 


(2.40) 

(2.41) 

(2.42) 


From  these  we  obtain 


S13 


1  +  cVu  I 
z 


(2.43) 


which  implies  that  the  effective  viscosity  for  simple  shear  is 
1 

2  2  2  * 

1  +  eVu* 
z 


(2.44) 
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The  next  steps  depend  on  which  of  the  constitutive  relations  (3=0 

or  3  “  1)  is  being  used* 

Upper  Convected  Maxwell  Model  (3=1) 

We  substitute  (3*6?)  into  the  constitutive  relation  (2*39)  to  obtain 

Cua  =  <Gh)(h-a)  +  X  (3.7P) 

where 

<*  a  *w  -  «"*  •  <>•»> 

Integrating  (3.7)  with  respect  to  z  and  UBing  the  boundary  condition 
(2,2(>)  we  obtain 

Cu  -  (Gh)(y2  +  ha  -  \  a2)  +  X  |j-  [  (Gh)(ha  a2)  }  .  (3.9P) 

We  Integrate  again  with  respect  to  a  over  0  to  h.  This  gives,  in  the 
notation  of  (3.2P), 

* 

ChQ  -  (Gh) ( y^h  +  |h3)  +  l(Gh)h3lt  +  ^h2[(Gh)h]t  .  (3.10P) 

We  now  substitute  this  into  (3,lP),  wh'ich  becomes 

ht +  ir  {ccbxA +  |h3>  +  n"  [<Gh>h3lt  +  lh2  [(Gh)hlt  \ m  0  *<3.up> 

This  is  to  be  solved  subject, to  the  initial  condition 

h(x,0)  -  hQ(x)  ,  (3.12P) 

the  boundary  conditions  (2*15P)  and  the  constraint  (2*17P).  In  addition 
we  aasume  that  h^(x)  and  h(x,t)  are  synunetric  about  x  -  0. 

Hie  axisymmetric  analog  of  (3.7P)  is 

OttB  -  (Dh)(h-a)  +  X  ~  {  (Dh)  (h  -  a)  )•  (3.7A) 

where 

Dh  -  ~  (hrr  +  7hr“Bh)  •  (3.8A) 
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After  two  integrations  and  application  of  the  boundary  condition  (2.26) 

we  obtain 

h(r,t) 

C  f  rudz-  r(Dh)  4  |  h3)  +fev  [(Dh)h3)t  +  £  rh2[(Dh)h]t 
0 

(3.10A) 

and  then  (3.1A)  becomes 

ht  +  r  4-  i-h3)  4-  t  [(Dh)h3]t  +  £  rh2  [(Dh)h]t}  -  0 

(3. 11 A) 

The  initial  condition  Is 

h(r,0)  -  hQ(r)  ,  (3.12A) 

and  we  have  the  boundary  condition  (2.15A)  and  the  constraint  (2. 17A) . 

We  assume  also  the  symmetry  condition 

dhn  ah 

d~  <o)  -  ft  <°*t>  "  0  •  <3-13A> 


Coro tat ional  Maxwell  Model  (P  ■  0) 

We  combine  (3.6p)  and  (2*43)  to  obtain 
Cu 


1  4-  eCu 


TTl  “  ^<sh)  <h  - z)  “  F  »  **y 


(3.14P) 


from  which  we  deduce  that 


Cu  - 

i 


1  -  A  -  4c2F2 
2e2F 


(3.15F) 


As  is  usual  in  flows  with  shear  thinning  a  solution  ceases  to  exist  if 
e  is  too  large*  Wc  assume  e  to  be  sufficiently  small  that  an  expansion 
of  (3.15)  is  possible,  whereupon  we  obtain 


Cuz  -  (Gh)  (h  ■  z)  4*  e2[(Gh)(h  -z)]3  4-  0(e4) 


(3.16P) 
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Integrating  and  using  the  boundary  condition  (2*26)  we  obtain 

Cu  -  (Gh)  jy2  +  j  h2  -  \  (h-z)2  }+  \  e2(Gh)3{  h4-  (h-a)4}  .  (3.17P) 

Integrating  again  over  0  to  h,  we  obtain 

ChQ  -  (GhKy2^  +  |  h3)  +  j  e2(Gh)V  .  (3.18P) 

Substituting  into  (3.1P)  we  now  obtain 

ht  +  ^{(GhXY^  +  ih3)  +  i  e2(Gh)3h5  }  -  0  .  (3.19P) 

For  the  axisymmetric  case  an  exactly  analogous  procedure  leads  to  the 
equation 

ht  +  7  b{r<Dh)<Y2h  +  3  h3)  +  5  e2r<Dh)3h5  }  “  0  (3*19A) 

with  conditions  (3.12A)  and  (3.13A). 

It  is  convenient  to  transform  to  a  moving  frame  in  which  the  position 
of  the  contact  line  is  fixed.  We  set 

5  -  x/a(t)  (3.20P) 

in  the  two-dimensional  case,  and 

5  -  r/a(t)  (3.20A) 

in  the  axisymmetric  case,  In  both  cases  the  interfacial  shape  is  given  by 

*-hfe,t)  (3.21) 

defined  on 

-1  <  5  <  1  (3.22P) 

or 

0  <  C  <  1  (3.22A) 

respectively. 
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For  the  motion  of  the  drop  the  equations  (2.18)  and  (2.20)  become 


2S-  (t)  -  C  [  -h?(l,t)/a(t)  -  eF]  (3.23) 

with 

a(0)  -  1  .  (3.24) 

The  various  forms  of  the  evolution  equation  are  transformed  as  follows; 
Equation  (3. IIP)  becomes 


£h  +  T  k  l  (6h)<v2h  +  J  h3)  +  12  fi[<Gh)h3]  +  £  h2£[(Gh)h]  }  «  0 

A 

(3.25P) 


where 


£  5  St  -  b<t>E  If 


b(t)  -  Ij-inaCt) 


Gh  ■  -  Ba2^  . 


(3.26) 

(3.27) 
(3.28P) 


Similarly  equation  (3.11A)  becomes 

£h  +  ~r  w  {Ufih>  <y2h  +  J  b3)  +  12  ^tCDbJh3!  +  |£-  h2Ct  (Dh)h] }  -  0 


where 


Dh 


3d  {  hd£  +  T  **d  "  B<l2b  |  * 


Equation  (3.19?)  for  the  corotational  model  becomes 

2 

*h  +  Tfr{ +  \ h3> +  H <6h>3h5 }  -  ° 

*  3a 


(3.25A) 


(3.28A) 


(3.29P) 


while  (3.19A)  becomes 

2 

Ch  +  |^-  {  5  CDbX-Y2!!  +“h3)  +S~d(Dh)V}  -  0 


ad 


(3. 29 A) 
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The  initial  conditions  are 


where 


b(5,0)  -  hQ(e) 


h0tt)  "  h0(x)  *  -1  <  x  <  1 


(3.30) 


(3.31P) 


h0(O  -  h0(r)  ,  0  <  r  <  1  (3.31A) 

and  the  constraints  (2.11) -(2.13)  apply  to  the  initial  shape.  The  boundary 


conditions  (2.15)  become 


h(-l,t)  -  h(l,t)  -  0 


(3.32P) 


h(l,t)  -  0  . 

The  volume  conservation  condition  (2.17)  becomes 


(3.32A) 


h(£,t)d£ 


(3.33P) 


in  the  two-dimensional  case,  and 


1 

-  2/eh(5,t)d5 

a  (O  Jn 


(3.33A) 


in  the  axisynmetric  case.,  We  Impose  also  symmetry  conditions  at  £  <*  0. 
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4.  SOLUTION  PROCEDURE 

The  problem  posed  in  the  previous  Section  is  soJved  by  transforming 
the  evolution  equation  (3.11)  or  (3.19)  into  a  system  of  ordinary  differ¬ 
ential  equations  in  the  time  t.  This  is  done  by  use  of  the  Galerkin  method. 
In  the  two-dimensional  case  we  define  functions 


h  (5) 
nsw 


(1  -5T  .  n  -0,1,2,... 


n  <  0 


for  integer  values  of  n.  Correspondingly  we  define 


(4.1) 


JL 

“»  '  / 


n  >  0 


(4.2P) 


These  are  then  given  by 


“0  ■  2  •  an  ■  ST+T  “n-l  •  "  L  1 


(4.3P) 


Next  we  introduce  the  Galerkin  expansion 
N+l 

h(e»t)  -  l  4>k<t)hlc(€) 


(4.4) 


where  the  ^  are  unknown  functions  of  t.  We  refer  to  (4.4)  as  an  N-term 
Galerkin  expansion,  since  one  of  the  $*8,  namely  Is  fixed  in  terras 

of  through  the  volume  conservation  condition  (3.33).  In  fact, 

substituting  (4.4)  into  (3.33P)  we  obtain 


Vw<‘> '  v/a<e>  -  J.W*)  • 

k«l 


(4.5P) 


All  the  boundary  conditions  and  constraints  on  h  are  satisfied  by  the 
representation  (4.4) -(4. 5). 

With  the  aid  of  (4.4)  we  construct  a  system  of  N  +  2  ordinary  differ* 
ential  equations  for  the  functions  <^(t) ,  •  •  •  »‘t,N+i(t)  ,a(t) •  First  we  have 
that  (3.23) -(3,24)  become 
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itoaiki.'  hjt V i 1  7, 


(4.6) 


with 


a  «  C  [  2<t^/a  -  9p] 


a(0)  »  1  . 


(4.7) 


Next  we  have  iron  differentiation  of  (4.5)  the  equation 


‘Wn+i  +  J3aA +  r " 0 


(4.8P) 


where  b  is  defined  by  (3.27).  Then  we  substitute  (4.4)  into  the  appropriate 
evolution  equation,  Multiply  by  hfl(5)  for  n  ■  1,2, • * • ,N  and  integrate  over 
(-1,4-1)  with  respect  to  g.  In  the  case  of  the  upper  convected  Maxwell 
model  evolution  equation  (3.25P)  this  gives  the  following  system  of  N 
ordinary  differential  equations: 

4<hn  3-{  (GhHySi  +  |  h3)  +  ~  C[(Gh)h3]  -r  £  h2£[(Gh)h]))-  0 


(4.9P) 


n  -  1..2,»**,N,  where  <■>  denotes  the  inner  product 


(Z)£tt)dK  . 


(4.10) 


We  now  have  N  4-  ?.  equations,  namely  (4.6),  (4.8P)  and  (4.9P);  they  are 
evidently  highly  nonlinear.  The  '’ondition  (2.12P)  on  the  initial  slope 
gives  the  initial  condition 


MO) 


(4.11) 


The  fact  that  we  wish  to  allow  for  the  drop  to  have  arbitrary  initial  shape 
is  expressed  through  the  initial  data 

4> 2(0) , •  •  •  j't’ jq(0)  arbitrary  .  (4.1,2) 
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Finally  we  have  from  (4.5p)  that 

N 

“N+l^N+l^0^  “  V  ‘  *  (4.13P) 

Thus  we  have  a  total  of  N  +  2  initial  conditions:  (4.7),  (4.11),  (4.12) 
and  (4.13P). 

In  the  case  of  the  corotational  Maxwell  model  (3.29P)  we  follow  the 
same  procedure  and  obtain  in  place  of  (4.9P)  the  system 

•4^'nf^+4ni{{(5,,)<Y2h  +  Jh3)+^6<6h)V}).0  ,  n  -  l,— 

5a 

(4.14P) 

The  other  equations  and  the  Initial  conditions  are  as  before. 

In  the  axi symmetric  case  we  again  use  the  system  of  basic  functions 
(4.1).  We  also  define 
1 

6n  "/»„«>«  •  »>°  (4.2A) 

ii 

which  gives 

_  1 

8u  "  2<fn  + 1)  »  n  -  0  *  (4.3A) 

The  Galerkin  expansion  (4.4)  is  again  used.  The  volume  constraint  (3,33) 
now  takes  the  form 

“N+lVl(t)  “  .  (4.5A) 

Equation  (4.b)  and  the  accompanying  condition  (4,7)  remain  applicable. 
In  place  of  (4.8P)  we  differentiate  (4.5A)  to  obtain 

.  n  •  Vb 

^W-l^Nfl  +  .1* k*k  +  ~2  "  0  *  (4.8A) 

k"l  Tta 
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The  evolution  equation  (3.25A),  alter  integration,  becomes 

a4  <Chn*h^  +  ^hn  ~  {  £(Dh)(Y2h  +  j  h3)  +  £[(Dh)h3J 

+  h2£[(Dh)h]|^  "0  ,  n  “!,•••, N  (4.9A) 

where  the  inner  product  is  defined  by  (4,10).  The  evolution  equation 
(3.29A)  becomes 

a'X'Vh)  +  <h„  If {  6(Dh) (Y2h  +  |  h3)  +  0  (Dh)Vj^> -  0 

n  ■  1,»**,N  .  (4.14A) 

The  Initial  conditions  (4,7),  (4.11)  and  (4.12)  continue  to  apply,  while 
(4.13P)  is  replaced  by 

<WW°>  -  h  -  j16k*k<°>  •  <413*> 

The  respective  systems  are  solved  numerically  using  a  modified 
standard  package  called  DGEAR. 
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